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Abstract 

We present a conjecture for the normalisation of the twist two conformal partial waves in a double 
OPE limit of the four-point function of stress tensor multiplets in M = 4 super Yang-Mills theory 
up to three loops. This contains information about the structure constants in the OPE. 

Like the twist two anomalous dimensions our result is expressed as a linear combination of 
harmonic sums whose argument is the spin of the exchanged operators. 

To arrive at the result we derive asymptotic expansions for the twist two part of two unknown 
three-loop integrals using the method of expansion by regions, complemented by some intuition 
gained on the example of the ladder integrals up to three loops. 



1 Introduction 

In M = 4 super Yang-Mills theory the loop corrections to the four-point function of stress-energy 
tensor multiplets take a factorised form: A superconformal invariant is multiplied by an x-space 
function expressed in terms of finite conformal integrals jU |2] . The one- and two-loop corrections 
to the x-space part have been known for a while [31 HJ E]. The three-loop contribution could 
possibly still be calculated using the methods of [4] based on superconformal invariance and 
M = 2 supergraphs, but such an endeavour would surely be rather cumbersome. 

However, it is not necessary to start from off-shell Feynman graphs in order to construct the 
loop corrections to the four-point correlator. One can rather sort the set of candidate scalar /-loop 
conformal integrals into orbits under an enlarged permutation symmetry discovered in [2], and 
then fix the overall coefficient of each orbit either by appealing to the correlator/amplitude duality 
[B] or by independent criteria relating to the suppression of the highest logarithmic singularities 
in accordance with the expected singular behaviour of the correlator in a Euclidean coincidence 
limit or in a light-cone limit, c.f. [7\. In the latter paper the integrand of the four-point correlator 
was constructed up to six loops relying on these criteria. 

In a Euclidean double OPE limit in which the positions of the four operators approach each 
other pairwise we can extract anomalous dimensions from a decomposition in terms of conformal 
partial waves [HJ. Here the two pairs of operators have an operator product expansion expressed 
as an infinite series of other operators, and the four-point function is essentially reduced to 
an infinite sum over the two-point function of these "exchanged" operators. In this paper we 
analyse the exchange of twist two operators. Their anomalous dimensions are given by a universal 
function depending on the spin. More precisely, perturbative calculations up to three loops have 
shown that this function is a linear combination of harmonic sums [HI UU] . This result has been of 
fundamental importance for the construction of an integrable system describing the higher-loop 
anomalous dimensions [11] in the so-called sl(2) or twist sector. 

The double OPE limit of the four-point function equally contains information about the 
structure constants in the OPE or equivalently the three-point functions (TTO^) where T is 
the stress tensor multiplet and is any (twist) operator occurring in the OPE0. We show in 
this article that up to three loops the structure constants of twist two operators are given by a 
universal function written in terms of harmonic sums. Recently there has been rising interest in 
an integrable systems explanation of structure constants [12]. We hope to foster this development 
with the formulae here presented. 

At one and two loops the quantum corrections to the four-point correlator are encoded in the 
functions 



9(1,2,3,4), 






+ 



2 (/i(l,2;3,4) + /i(l,3;2,4) + h{l,4;2, 



3)). 



This picture does not require an explicit definition of the quantum corrected operators. 
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At three-loop level we use the result from [2] 

= 2 0(1, 2, 3, 4) (x 2 12 x 2 34 h(l, 2; 3, 4) + x 2 3 x 2 M h(l, 3; 2, 4) + x 2 u x 2 3 h(l, 4; 2, 3)) (3) 

+ 6/(l,2;3,4) + 6/(l,3;2,4) + 6/(l,4;2,3) + 4 J E(l;3,4;2) + 4 J E(l;2,4;3) + 4 J E(l;2,3;4) 

+ (x 2 3 x 2 24 + x 2 *x\ 3 ) H(l, 2; 3, 4) + (xj 2 x 2 u + x 2 u x 2 23 ) if (1, 3; 2, 4) + 
+ (x 2 12 x 2 u + x 2 13 x 2 2i )H(l,A;2,3). 

This is the full answer; there are no non-planar corrections up to this order. The definition of 
the integrals is 

(1 2 3 4) = — [ 

ft7I J X 15 X 25 X 35 X 45 
hf-t n.o a\ 2 1 /" <j 4 X 5 G? 4 X 6 

rt^l, Z, 6, 4) - X 34 , . 2 / 2 2 2 x 2 f~2 J J \ ' W 
^71 J J \ X 15 X 35 X 45/ ,X 56\ X 26 X 36" ( '46J 



/(1,2;3,4) = x 



1 /" d x 5 d x e d x 7 

/I tt-2 ^ 3 / ( ,->-» 2 ™ 2 ™ 2 \ ,>.2 / , 2 ™ 2 \ ,y»2 / ™ 2 ™ 2 ™ 2 

^71 ; j yx 15 x 35 x i5 )x 57 yx 37 x 47 )x 67 yx 26 x 36 x t 



£7(1, 2; 3, 4) — x 23 x 2i r^zyi 



m) 

d A x 5 d 4 x 6 d A x 7 x\ 7 



( /1'tt-2A3 / ( />^2 ,->^2 ry»2 \ ^y»2 / ,>^2 ,-y»2 ,-y»2 \ ,-y»2 / ,-y»2 ,->•» 2 \ 
ytn ) J \X lh X 25 X 3 e > )X§ 7 \X2 7 X 31 Xtf)X$ 7 \X 1 QA 2 Q- L 46) 

1 f d x$ d A XQ d A x 7 x\ & 

(,-fi 2 ,-y-»2 ,-y,2 ™2 \ ™2 / ™2 ,-y,2 \ ,-y»2 / ^2 ™2 ,-^,2 ™2 ^ 
X 15 X 25 X 35 X 45J X 57\ X 37 X 47J X 67\ X 16 X 26 X 36 X 46J 



In (|3]) we encounter two conformal non-ladder three-loop integrals E,H (for "easy" and "hard") 
which are not explicitly known. Fortunately, one can extract the twist two trajectory directly 
from the integrals by the method of "expansion by regions" [13]. To this end we move one 
point to infinity — which is always possible due to conformal invariance — and then study an 
expansion of the resulting three-loop integrals in terms of the remaining small distance of the 
OPE limit. The procedure collapses the scalar three-point integrals to two-point tensor integrals. 
We evaluate the genuine three-loop pieces by tensor reduction and an integration routine of the 
Mincer system [T4]. The information so obtained is not quite sufficient for our purpose so that 
we have to supplement it by some intuition about the form of the results. 

It is in principle possible to extend our analysis to higher twist but we would meet more and 
more difficulties in extracting enough information from the integrals using Mincer. 

The article is structured as follows: In Section 1 we derive an asymptotic expansion of the box 
integrals in terms of harmonic sums with only positive indices. Section 2 discusses the method 
of expansion by regions and our results for the twist two part of the unknown integrals E, if. In 
Section 4 we discuss the conformal partial wave decomposition of the correlator. 



2 Asymptotic expansion of the box-integrals 

The massless L-loop boxes (equivalently the (L + l)-rung four-point ladder diagrams) have been 
evaluated in [15]. They are conformal integrals given by the appropriate weight factor times 
$( L )(x, y) with the functions 

L-l 



log*" 1 ® (log (I) + log(O) ~ (log (!) + 2 log(0 
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The arguments are space time cross ratios. In a Euclidean coincidence limit x\i, X34 — > 



22 22 
, = 1-7 = ^^1, M= V|4^ . (6) 

X 13 X 24 J '13 x 24 

and similar with 23 <H- x 4 . Due to the x -H- ?/ reflection symmetry of we may analyse all 
four such situations by considering the limit x — y 1, y — > of Since we will be interested in 
the y° part of the expansion we can conveniently put y = in the first factor (the denominator 
under d£) of the integrand; the parameter integrals remain well-defined. The simplest expansions 
arise for g(l,4,2,3), /i(l,4;2,3), /(1,4;2,3) where we put x = 1/v, y = u/v. The ^-integration 
leads to Li n (—Y) which is replaced by its Taylor series for small Y. To take the same limit on 
/i(l,3;2,4), Z(l,3;2, 4) we choose x = v, y = u. The integration then yields Li n (—Y/(1 — Y)) 
which can also straightforwardly be expanded in Y. 

The point exchange 23 <H- £4 implies the transformation 

u , Y -> — (7) 

i-y' i-Y v ; 

on the cross ratios. Combined with a division by 1 — Y due to the outer weight factor this 
exchanges the expansions of h(l, 4; 2, 3), 1(1, 4; 2, 3) in the Xi 2 ,x 34 — >■ limit with those of 
/i(l,3;2,4), Z(l,3;2, 4), respectively. The one-loop box #1234 is totally symmetric, and indeed 
the two choices for x, y yield the same series. 

The third case, g(l, 2, 3, 4), h(l, 2; 3, 4), 1(1, 2; 3, 4) corresponding to the choice 

x = - , y = - (8) 
u u 

(or its reflection) is mapped into itself. It cannot be analysed by the same simple manipulation 
of the integrand in (E]). For the fully symmetric one- loop box this poses no problem because the 
series expansion must fall upon what we had before. In order to analyse this limit of the higher 
box integrals we write the first factor of the integrand as 

u ( d£ d( \ ^ 2-Y-u±\ x — -r 

A U Y i \ ) ' fa = 2(l-y) ' *=V-4«+ (!- + «)'. W 



Using 



we find 



-L = l-z±, x± = ^(F + n±A) (10) 



- $ (2) (y, u) = -6 (Lz 4 (l - x+) - Lu(l - xJ)) + 3 log(l - Y) (Li 3 (l - x + ) - Li 3 (l - 



u 



-~ log 2 (l - F) (Li 2 (l - x + ) - Li 2 (l - x_)) , (11) 
A $ (3) (F, it) = -20 (Li 6 (l - x+) - Li 6 (l - + 10 log(l - y) (Li 5 (l - z + ) - Li 5 (l - x_)) 



-2 log 2 (l - y) (Li 4 (l - x+) - Lu(l - z_)) 

+1 log 3 (l - Y) (Li 3 (l - x+) - Li 3 (l - x-)) . (12) 
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We seek an asymptotic expansion in a Euclidean regime where Y,u > 0. In the coincidence limit, 
Y 2 and u are of the same order in small quantities. Elementary trigonometry shows —Au+Y 2 < 
to leading order, so that the root A is purely imaginary. Hence x± are complex with $l(x±) > 0. 
After expanding Li n (l — x±) in x± in this regime we treat the variables Y and u as independent 
and further expand first in Y and then in u. The individual polylogarithm terms contribute 
Y n /u m at order n >= 0, with m = [n/2] + 1/2, [n/2] , ... up to positive powers of u, but negative 
and half-integer powers of u cancel in the complete expressions for <f>( L \ We retain only 0(u°), 
the twist two trajectory. 

When the branch point is approached from the left the polylogarithms are described by the 
following asymptotic series: 



Li 2 (l — x) - 
Li 3 (l — x) - 
L« 4 (l — x) - 
Li 5 (l — x) - 
Li 6 (l — x) - 
with the functions 



c+E- 

\ - x n 

1 

oo 

1 

oo 

C5+E 

i 

oo 

<»+E 



X' 

n 

x n 
n 

x n 
n 



x 
n 



- - \og{x) Iq + Ji 



<2lo 



- - \og{x) h + J 2 
n 



I (.2 h ( - - l0g(x) ) I 2 + J 3 



C4/o-C3/l-C2/2- 

-C5 lo ~ U h ~ C3 h ~ C2 h - I - - lo g(^ 



log(x) ) J 3 + Ja 

1 



(13) 



and 



/o = 1, 
/1 

h 



h 
h 



-^2 + ^! 2 , 



2 
1 



3 2 6 1 

1 1 1 2 1 2 1 4 

— — S4 + — Si S3 + — 5*2 ~~ t S1S2 + — S^ 
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24 



(14) 



Ji = 0, 

J2 = S2 , 

^3 = 5*3 — S1S2 , 



J4 — 5*4 



Si S3 



2 
5 



Si 



S\ 5*2 , 



(15) 



1 2 1 2 1 3 
^5 — S 5 — S1S4 — - S2S3 + - S X S3 + - SiS 2 — - S 1 S 2 . 

o 2 I o 

In these formulae the range of the harmonic sums is from 1 to n, so they denote Si(n), S 2 (n) etc. 
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The expansions are of pure transcendentality: A sum ^2 in x n /n m = Li m (x) is obviously 
assigned weight m. The more general expression Y2 n x n /n m Si(n) will be regarded as a weight 
I + m object if Si is any harmonic sum (or a product thereof) of total weight I. The redundant 
symbols I Q and J\ were introduced only to emphasize the pretty iterative pattern, by which the 
coefficients of (^-values in the asymptotic expansion of Li n (l — x) is known from the expansions 
of Li m (l — x) with m < n. The I n , J n functions that we display here have been matched on 
Mathematica output. Unfortunately, Mathematica runs into problems when the order of the 
expansion becomes too high. The iteration as well as the fact that I n , J n can apparently always 
be expressed as products of simple values were observed on the lower examples in the list, 
where the fits can contain only a very limited number of distinct structures. Once a fit has been 
established it is a trivial matter to continue the original series up to very high orders. 

The asymptotic expansions for the entire box integrals (the complete functions including 
rational pre-factors) are all similar to (fT3"j) . We have begun by analysing the highest logarithms 
and/or Q- value contributions to understand what type of object to fit to the series. 

By the simple manipulation on the integrand of (J5J) skeched in the beginning of this section 
we found at 0(u°): 



4 ,7(1, 4, 2, 3) -> 



4^(1,45 2 , 3) 



n=l 



n 



4^(1, 4; 2, 3) -> J2 



Y' 



n=l 



n 



E 

n=l 



Y 



n—1 



n 



log(w) + - 
n 



1 i 2/ \ 1 i / \ 3 6 
- log [u) log(w)— + — 

2 n n z n 6 



10 20 



lo S ( u ) -o + o lo S ( u ) ~1 - lo g( M ) ~i + 
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rr 



n 4 n 5 



(16) 
(17) 



The limits of 77.(1, 3; 2, 4), 1(1, 3; 2, 4) are of a more complicated form: The resulting series can be 
fitted on linear combinations of harmonic sums with exclusively positive indices: 



4 h(l,3; 2, 4) -> 



Y 



n-l r 



n 



E 

71=1 

4Z(1,3;2,4) -> 



log 2 (u) 5*1 - log(u) 



Si 



n 



2S- 



E 

71=1 



Y 



n-l 



II 



log(u) 



S 2 

2n 2 ^ 2n 2 



?1 + ^ _ + 2 S 3 + 2 Si 2 



77" 



77 



S 2 S2 

- — I- - — h 5*1 02 + S3 + Si 2 
2n 2n 



q CO o o 

02 , Oi 02 03 . ^1,2 a 2 



+ 



77 



n 



n 



+ 5*2 + S\ S3 + 2 S4 + 2 Si^ + Si, 1,2 — Si^, 



(19) 



(20) 



02 C C C C C C2 C C O C O C C C 

J I <->2 Si 02 S3 »->l,2 1 £2 1 ^L_i 1 Z_£ 1 _ X ' 3 1 ^l.l.S _ 13 1,2,1 



2n s 



2n s 



77^ 



n- 



77- 



77 



77 



77 



77 



77 



77 



+ 2 S2 S3 + Si 04 + 3 S5 + 3 Si^ + 5*2,3 + 2 Si, 1,3 — 2 Si i3j i + S2,i,2 — S2,2,i — Sx,i t 2,i + 512 .1, 



In these equations — like anywhere in this article 
explicitly stated otherwise. 



all harmonic sums have argument n unless 
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The coincidence limit on the remainig cases h(l, 2; 3, 4), 1(1, 2; 3, 4) requires the more compli- 
cated procedure outlined above. We can match the results at O(u ) by the following expressions: 



x 



13 



fc(l,2;3,4) -> 

n+1 



n=0 

oo 



x 



13 



n=0 



6C(3) + 5 li2 - 
20C(5) + C(3) (^ 2 -5 2 ) 



(21) 



"32,3 + £3,2 + Ox, 1,3 — 5*3,1,1 — 5*1,2,2 + 5*2,2,1 — 5*1,1,2,1 + 5*1,2,1,1 



Although this is not manifest, these expansions — as well as the simple result for g(l, 4, 2, 3) - 
are mapped onto themselves under the point exchange x 3 O x 4 . 

3 Limits of E and H by Asymptotic Expansion 

The non-ladder integrals in the three-loop correction to the four-point function are not yet 
explicitly known. Fortunately, the method of "asymptotic expansion of Feynman integrals" 
[13] allows us to analyse Feynman diagrams in any limit; coincidence limits on finite Euclidean 
integrals are almost the defining examples. Like in the case of the box integrals we will obtain 
the leading terms of a power series in Y at u° and seek a fit on harmonic sums. 
Recall the definition 



£(1,3; 2, 4) 



d X§d Xq(1 Xj ^23^34^17 
2 2 2 2 2 2 2 2 2 2 2~~ 

iy" nn" /y nn" nn" nn" /y /y" /y" iy" nn" 

x 15 x 16 x 25 x 27 x 46 x 47 x 35 x 36 x 37 x 57 x 67 



1 



/-V>2 ,-y»2 

x 13 x 24 



U,V) 



(22) 



where the second equality follows by conformal covariance. The integral representation is in- 
variant under the exchange of points 2 and 4. As this exchanges u and v we conclude that 
& E \u,v) = & E \v,u). Next, the exchange of points 1 and 3 has the same effect on the cross 
ratios, whereby we can conclude that £'(1,3; 2, 4) = £(3,1; 2, 4) despite of the apparent asym- 
metry of the integrand between x\ and X3. All in all we are left with three cases to analyse: 
£(1,2; 3, 4) and £(1,3; 2, 4), £(1,4; 2, 3). In the limit £12,234 — > 0, the first case stays apart, 
while the other two are related by the exchange of points 3 and 4. 



3.1 E(l, 4; 2, 3) and £7(1, 3; 2, 4) 

Due to conformal covariance the integrals can be uniquely reconstructed from a limit where, say, 
point 4 is moved to infinity by replacing 



rf ry 2 I V / > ■ 2 

x 12 x 34 " x 12 1 



2 2 2 
x 13 x 24 ^ x 13 J 



iy ^ iy ^ J k iy *^ 

x 14 x 23 ^ x 23 



Note that in the limit 



u -> 



X 



X 



12 
2 ' 
13 



X 



23 



2(xi 2 .xi 3 ) 

1 K h U . 



X 



13 



X 



13 



We consider 



£7(1; 2, 3) = lim x\ £(1,4; 2, 3) = / 

X4— >oo J 



d x^d x^d^Xj x 2 



17 



22222222 

ry' fy** ty*- 1 nr*** np" iy** rf" 

x 15 x 16 x 25 x 27 x 36 x 37 x 57 x 67 



(23) 



(24) 



(25) 
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Thanks to translation invariance we may further put X\ = 0. Let us re-label x 2 — Pi, x 3 = 
p 2 , x 5 = ki, Xq = k 2 , x 7 = k 3 . The integral becomes 

f d A kid A k 2 d A k 3 kj 

{Pl,P2 ' = J k\k\{k x - k 3 )\k 2 - hfih - Pl f{k 3 ~ Pi) 2 {k 2 - p 2 f{k 3 ~ P2) 2 (2b) 

and we are interested in the limit p\ — > 0. We might simply try to expand the integrand in pi 
using 

1 = if ( 2(k 1 . Pl )-pi yi 1 = 1 f> ( 2(h. Pl )-pi y 2 
ih-viY k\^\ k\ J ' {h- Pl y ki^ o { ki ) ■ w 

These equations are only valid if k\,k 3 > p\, of course, and for example in calculations with 
orthogonal polynomials one would indeed subdivide the integration domains according to the 
validity of such expansions. Here we rather put the measure into D = 4 — 2 e dimensions 
in order to regularise the IR singularities l/(k\) mi and l/{k 3 ) m2 that arise by extending the 
integration domain to the origin. We have expanded the integrand according to the formal 
assigment ki,k 2 ,k 3 = 0{p 2 ) » p\. We call this the "top region". In the "bottom region" we 
declare ki,k 2 ,k 3 = 0(pi) « p 2 and likewise employ the geometric series to expand the k 2 — p 2 
and k 3 — p 2 propagators in k 2 ,k 3 , respectively. Here we find UV poles arising from the part of 
the integration domain where the momenta are large. The method of "expansion by regions" 
consists of evaluating not only the top and bottom regions, but the sum of all eight possibilities 
arising from ki = 0{p\) or 0(p 2 ). All singularities cancel and the logarithms combine into powers 
of logw. The expansion by regions is equivalent to the "expansion by subgraphs" which in turn 
has been proven by renormalisation theory to yield valid asymptotic expansions [T5] . 

In massless theories, in some regions one encounters "no-scale" integrals J d D k/(k 2 ) a = 0. 
In the case at hand we find non-vanishing contributions in the regions 

Rx : h,k 2 ,k 3 ~ p 2 , R2 ■ k 2 ,k 3 ~p 2 - k x ~ Pl , (28) 

#3 : h ~ p 2 ; fei, k 3 ~ p x , R 4 : fci, k 2 , k 3 ~ p x . 

Here R\, R4 are the "top" and "bottom" problems mentioned in the last paragraph. Quite 
generally, the top problem is an /-loop integral with very high indices (exponents) on some 
propagators and therefore the hardest to solve. The bottom problem is often of the same topology 
but the indices are all as in the original integral. The mixed cases break into an m-loop integral 
depending on the small scale and an (/ — m)-loop integral depending on the large scale. 

The original integral E{1\ 2, 3) has dimension [1/p 2 ]. Let us consider the bottom problem R4. 
To lowest order we put k 2 —p 2 ,k 3 —p 2 — > p 2 . It follows that the corresponding three- loop integral 
must go as p\. We conclude that the bottom problem does not contribute at order u°. The same 
applies to R 3 : To lowest order we have l/(k 3 — p 2 ) 2 —> l/p\. We find a one-loop integral that 
has dimension of l/(pl) 1+t , so once again the two- loop integral depending on pi must produce 
at least one power of u. 

The leading term of the top problem has p- dependence ^-/{p 2 ) l+3e , thatofi^ is l/((p|) 1+2e (Pi) e ) 
Since there are only two contributing regions we should find that the leading poles from Ri, R 2 
are equal and of opposite sign. Further, a pole e~ n yields logarithms to the n-th power in the 
finite part of the e expansion. The logarithms from the two regions can only combine into log(w), 
if the leading singularity is a simple pole. 
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We are finally in a position to consider the top region Ri in detail. We can drop p\ from ([2 
because we want to restrict to 0(u°). We wish to calculate 

F( ,to P r d D hd D k 2 d D k 3 (2(k l .p 1 )r^2(k 3 . Pl ))^ 

[PUP2) U - 2^ J ( A; 2)2 + n l(A . 1 _ fe )2^^2 ) „ 2 ^ 2 _ fe)2(A . 2 _ p2)2(A . 3 _ p2)2 - I ) 

In terms of the classification of three-loop propagator type integrals in [TJ] this integral is of 
topology 2 , the two- loop master 7\ with a bubble insertion in an outer line. The k\ integration 
(i.e. the momentum going around the bubble) can be executed using 

f d D k P s (p) __*__xrn( R ■J(^)V D *V pml ( ^ 

J wrW^pff = (fr^ g^'^prUJ m) ) k _ p m 

where P s (k) is any polynomial of order s and = d 2 / dk^dk^. Last, 

r(a + /3 + z-D/2)r(D/2-a + g -Qr(D/2-/3 + 

G(ot,P,s,z) = — — . 31 

r(a) r(p) T(D — a — p + s) 

In our case the numerator polynomial is k^ 1 . . . k^ 711 , all contracted onto p\. In applying (130|) to 
( 129]) we discard terms with i > since they yield p\. It follows 

uod „, ^ /" d D k 2 d D k 3 (2(k 3 .p 1 )) nM . . 

The region gives a very similar result: We expand 1/(^3 — £>i) 2 in pi and 1/(^3 — &i) 2 in A^. 
Powers of k\jk\ can be dropped because they vanish under the k\ integral, which is a one- loop 
bubble with ingoing momentum p\. It multiplies a T\ topology depending on p 2 . Using (|30|) for 
the bubble we find 

d D k 2 d D k 3 (2(A; 3 .pi)) ni+n2 



1 r 



k 2 (A; 2)i + n 1+ n 2(A . 2 _ fc 3 )2(fc 2 _ p 2 )2(fc 3 _ ^ 



(33) 



n.i,n2=0 

From the last two equations it is already clear that the leading poles will cancel: 

(7(2 + 711,1,711,0)1^ = = -(7(1,1,7^,0)1,-! (34) 

e n + 1 

Recall that the first two arguments of the G function label the indices of the propagators in a 
one-loop bubble integral, reflecting the fact that (7(2 + m, 1, Tii, 0) has a pole of IR origin while 
(7(l,l,ni,0) contains a UV divergence. This is a very direct illustration of the cancellation of 
singularities due to extension of the integration domains to the whole of Minkowski space. The 
remaining two loop integrals in (132]) and f )33|) are equal to leading order in e. The traceless part 
of the numerator (i.e. no pi) always leads to finite integrals so that the recombination of the 
logarithms can take place as we had anticipated. 

Summing up, conformal invariance has permitted to reduce a four-point integral to a (generic) 
three-point one. The strategy of "expansion by regions" makes it possible to extract a power 
series in a small ingoing "momentum" p\ = x 2 i whereby the three-point integral collapses to a 
collection of two-point problems on the expense of introducing numerators with open indices. 

Even if there are numerators, the "rule of the triangle" [H] can be used to solve the two- loop 
master topology T\ with all integer exponents as in ( 1331) or one non-integer exponent on an outer 
line as in fl32|) . The triangle rule is an IBP identity 2 ] which reduces the exponents of one of 



integration by parts 
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the &2, k2 — k^, k% — P2 lines in favour of increasing those of k^, k^ — P2- In the case at hand a 
single application of the triangle rule would be sufficient because the lines involving fc 2 all have 
exponent one. Hence one of them will be cancelled in every term upon which the integral can be 
evaluated by iterated use of equation (|3"U|) . 

In practice we rather implement a different strategy: We will use tensor reduction in order to 
pull pi off the integrals and then apply the rule of the triangle to the resulting scalar integrals. 
In the top and bottom problems of the expansion by regions in other limits on E and H we 
encounter the three- loop FA topology which is implemented in the powerful "Mincer" system 
[H]. Since this programme only deals with scalar integrals we have to understand the tensor 
reduction at any rate. 

Let an Z-loop integral depending on a single outer scale P2 have a numerator jV^ 1 "^ 8 with s 
open indices and some denominator D. Upon integration 



d u k 1 ...d u k l N^- 



Po Ms (p 2 )/o(p 2 ) + ... + P^/ s (p 2 ) (35) 



jj — 1 \f2J m\F2l [ s /2] 



Pf ^(Pi) = 4^(33.) (30) 



where 



As in ( 130]) the d'Alembertian replaces p^P^ D Y in a totally symmetric fashion. Symbolic 
differentiation gives a fast and simple algorithmic realisation of the tensor decomposition because 
we need not explicitly symmetrise. 

In order to determine the U in (l35j) we project the entire equation with P (p 2 ), . . . , P[ s /2](p2) 



obtaining 



■/ = M lt Ij . M,j = { ^ (5*Y Ml (3* ) U,i,r } . (37) 



4 



=P2 



Jt j m (a k y j- d D k 1 ...d D k l k" ...k" .\ ; ,,.,, , (gf 

The matrix My depends on the regulator because d p .p = D = 4 — 2e. However, the tensor 
reduction procedure is well-defined also in exactly four dimensions, so that M must have a power 
series expansion in e. To construct the inverse up to and including e 3 we only have to invert the 
lowest term and complete it by linear perturbation theory; in practice the whole procedure was 
still fast up to spin 80 in the most straightforward implementation. In addition, we can limit our 
scope to the construction of I because the trace terms in Piip?), i > obviously yield powers 
of u upon multiplication with p\. Finally, we point out that all steps of the algorithm can be 
performed on the numerator N before calling any integration routine. 

In our results we can finally replace pi,p2 by £21, £31, respectively, and then recover the 
dependence on the fourth point by identifying 2(xi2.xi3)/x^ 3 = Y + 0(u) and extending the 
overall denominator to x\ z x\ Al c.f. ( 123]) and f T24l . The limit of E(l, 3; 2, 4) can be obtained from 
this expansion by the transformation ([7]) corresponding to the exchange of points 3 and 4 and 
division by (1 — Y) owing to the overall denominator. The asymptotic series so obtained have 
very much the same features as the limits of the box integrals. By a short Mathematica script 
for the evaluation of the T\ integrals via the triangle rule we could easily generate the asymptotic 
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expansion up to order Y 31 , which is sufficient to determine all coefficients in fits on harmonic 
sums with only positive indices (possibly divided by powers of their arguments): 



4, £(1,4; 2, 3) 



(39) 



E 



Y 



n-l 



n=l 



n 



tog(ti)C(3)(6 5 1 )-tog(tt)f^-- 



Sf S\ S 2 



n 



<J1,3 + <33,1 + 2 Si,i |2 — 2 5*2,1, 



+ C(3) ( —-4^+12^ 



2 Sf 2 Si 5*2 2 So 2 Si S3 



n 



4, £(1,3; 2, 4) 

Y"n— 1 



- 4 S M + 4 S 4 ,i + 4 Si,i, 3 - 4 5 3 ,i,i + 4 S li2)2 - 4 5 2 , 2 ,i 

(40) 



E 

n=l 



n 



\og{u) C(3) - - log(« 



3S 4 



+ 3 S 1)3 - 5 1)2:1 + S 2jl)1 + C(3) + — - 2 S 2 



2S X 


5? 


n 3 


2n 2 




45 




+ 


V n" 1 


n 



So 



2 5*3 Si, 2 Sf 

1 Oi D 3 

n n 2 
8 Si 5 Sf 3 S 2 



n' 



2n 3 2n 3 



Si S2 3 S3 Si, 2 Sf Si S3 8 S4 4 Si, 3 Si, 1,2 _ 3 Si, 2,1 2 S 2 ,i,i 
n 2 n 2 n 2 n n n n n n n 

— 5 Si S 4 — 4 S2 S3 — 9 S5 + 13 Si,4 + 5 S2,3 — 2 Si, 1,3 — 2 Si,3,i + 4 S3,i,i 

+ "2,1,2 — S 2 ,2,l + "1,1,2,1 ~~ "1,2,1,1 

3.2 #(1,2; 3, 4) 

This is a single-log limit of the H integral. The derivation of an asymptotic expansion for the 
twist two trajectory is very similar to the single- log limit of the E integral, just that here the 
23 -H- 24 exchange maps the expansion to itself. In the limit £4 — > 00 we obtain the reduced 
integral 

#(1,2; 3) = lim x\ #(1,2; 3, 4)= / 

X4 ^°° J A 15 J '16- X '25- X '26 J '35- X '36- X '37- X '57 A 67 

We put £5 = k\, xq = k 2 , xj = fc 3 , x\ = 0, X2 = Pi, £ 3 = ^2 as before and expand using the 
geometric series according to which "loop momenta" are supposed to be of order p\ (small) 
or order p 2 (large). There are two contributing regions: The top problem ki ~ p 2 and k\ ~ 
Pi, k 2 ,k 3 ~ p 2 (and symmetrically k\ -H- fc 2 ). The bottom problem is non-vanishing but only 
comes in at 0(u). 

The top region yields an FA topology, with high indices for two propagators (p4,ps w.r.t. 
the definitions of [H]). The second region yields a one-loop bubble times the two-loop master 
Ti quite as R 2 in the evaluation of E in the last section. The Mincer system could derive the 
expansion of the top problem only up to and including 0(Y 20 ). This may point to an installation 
problem, but more likely this means that the recursion for the evaluation of the T\ master with a 
non-integer exponent on the central line (which is internally encountered upon applying the rule 
of the triangle to the FA topology) is not tabulated to sufficiently high orders. We obtain an 
asymptotic series containing \og(u) C(3), log(w), £(3) and purely rational terms. The coefficient 



d^Xscftxed^x? x\ % 

~~2 2 2 2 2 2 2 2 2~~ 

,-yiij ry-> np^ /y^ ry>^ /y>" nr>" nn" rr>^ 



(41) 
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of log(u) C(3) must have transcendentally weight two, and from the numerators we could in fact 
easily recognise Y n ~ 1 Si(n) /(n+1). Fits of this type are quickly found also for the non-£ log term 
and the £(3) term without logarithm. However, the expansion up to spin 20 does not furnish 
enough information to fix all 32 constants in our ansatz Si(n)/(n + l) 5 , . . . , S^(n)/(n + 1), . . . for 
the purely rational part (all sums were taken to have positive indices only). 

To make progress we invoke invariance under x 3 -H- x\ acting on the expansion 

oo 

lim H(l, 2; 3) = V Y n ~ x (log(u) a n + b n ) + 0{u) . (42) 

n=l 

On the cross ratios we have the aforementioned transformation ([7]). It remains to divide by 

(1 — Y) 2 , to re-expand in Y and to equate with the original series. The resulting system of 

equations allows to express one half of the a n , b n in terms of the others. We substitute the ansatz 

SAn) . Si(n) , S 5 (n) J _, , Soin) , , 

a n = c 5 ,i-^ + ... + C2,iC(3)^7t, b n = ^^ + ...+03,1( 3 ^r + ... 43 
n+1 n + 1 n + 1 n+1 

into the first, say, 60 such equations finding seven linear relations between the 15 constants 05^, 
one for three 03^ and fifteen equations on the set of 32 constants c^. The actual power series 
derived from the integral is consistent with these conditions and now suffices to pin down the fit 
at 0(u°): 

o4 #(1,2; 3, 4) -> (44) 



y — 

^ n + 1 

n=l 



- log(u) C(3) (24 5i) - log(«) (-2 S 2 2 + 4 S 1 S 3 + 2 S 4 - 4 S h3 + 4 5 1)1>2 - 4 



+ C(3) [^-6S! + 6S 2 )-^- + + _ Ml£ + 8 ^ 8 S W 



^n + 1 / n + 1 n + 1 n + 1 n + 1 n+1 n 

2 5*2 ^3 + 8 5*1 £4 + 10 S5 — 8 Si t 4 — 12 6*2,3 + 10 — 8 Si,3,i — 2 53^^ 

2 Si ,2,2 + 2 02,2,1 — 2 Si 1,2,1 + 2 SV 2 

,1.1 



(45) 



3.3 £7(1,2; 3, 4) 



We send point 3 to infinity and identify xi = pi, x 2 = 0, x 4 = p 2 , x 5 = ki, x 6 = k 2 , x 7 = fc 3 . 
The integral to expand is 

F( ]= 2 f d D kid D k 2 d D h (kj-2k 3 .pi+p 2 1 ) 

{Pi,P2) P 2 J fc fA;iA ; I(A ;i -A ; 3) 2 (A ;2 -A ;3 ) 2 (A ;i -p 1 ) 2 (^-Pi) 2 ^i-P2) 2 (^-P2) 2 1 J 

There are six contributing regions: The top and bottom problems, two regions with two integra- 
tion momenta treated as large and one as small, and two for the opposite. One of the latter does 
not contribute at O(u ). The top problem is of topology 2 , the bottom one of topology FA and 
the three other regions give (nested) bubble integrals. For the top problem and the trivial regions 
we need not appeal to Mincer, for the bottom problem we used its inbuilt FA routine once again. 
Unfortunately the programme was not able to go beyond Y 17 . To make matters worse, the log 3 
part of the asymptotic expansion is easy to guess and one sees that there are series of the type 
y n_1 S 2 (n)/n and y n ~ 1 S 2 (n)/(n + 1), whereby there are twice as many constants as before. 
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On the other hand, the expansion by regions must make sense also if the three terms of the 
numerator are treated separately: E(p 1 ,p 2 ) could be a generic three-point integral, for which 
the strategy should be operational independently of the fact that it is a limit of a conformal 
four-point integral. Interestingly we find that the k\ numerator is linked to the 5(n)/(n + l) type 
expansion while the other pieces cause the S(n)/n part. Consequently, running the expansion 
by regions separately for the two parts of the numerator we can double the number of conditions 
available from the integral itself. 

The integral must have point exchange symmetry x\ — > x 2 . This maps the 5(n)/n part of 
the fit to both pieces, but fortunately the S(n)/(n + 1) piece is sent onto itself. In this sector 
we could solve up to five constants. Next, in all the combinations of harmonic sums that we 
have encountered up to now the coefficients are integer or half integer. Using this knowledge one 
can simply play through all values for the remaining undetermined coefficients in a typical range 
and see where the dependent constants come out with denominator 1 or 2. Admittedly, this is a 
somewhat experimental approach, but there is very clearly only one reasonable solution: 



4, £7(1, 2; 3, 4)-> 



(47) 



E 



Y 



n-l 



n=l 



n 



■g log H 



+ - io g 2 H 



2Si 
n 

i2 



si 



4Si 3S 1 3S 2 „ „ „ 

H ~ 0± 2 — O Ol 2 



rr 



2n 



2n 



log(«) 



6Si 2Sl AS2 3 5i5 2 5s 

o o 



2SV 



3 5*9 5 5 4 
— ~~ + — — 



n J n* n n n 

5 51,3 - 51,1,2 + 251,2,1 - 5 2 ,i,^ + 20 C(5) + C(3) 
1 •.> c - c c •> c g i2 



2 Si 



n 



-Sf-5 S< 



85i_55f 13 5 2 5 Si S 2 2 S 3 



2 Si 5 3 5 S? 



2n s 



n* 

2 5i,i,2 3 5l,; 



2n A 



n" 



n 



3 Sa 

— - H 

2n 2n 



5 5i, : 



n 



_SiXi_ s 2 S 3 + 9S 5 -AS 1A -AS.. 



2,3 



n n n 

"5i,i,3 + 4 5i,3,i — 3 53,1,1 + 5i,2,2 — 2 5 2 ,i, 2 + 52, 2,1 + <Sl,l,2,l — 5i,2,l,l 



+ 



E 



Y 



n-l 



n=l 



n + 1 



-ilog 3 (5 1 2 + 3 5 2 )+ilog 2 ( M ) 



2 6 So , _, , _, 

+ —^7+453 + 45, 2 



log(u) 



45 2 



+ 



12 5, 



8 5 3 
+ -: + 



n + 1 n+1 
8 Si 9 



(n + 1) 2 (n + 1) 2 n+1 n+1 



- 3 Si 5 3 + 5 2 2 + 2 5 4 + 8 5i, 3 



+ 5!,i, 2 - 2 5 1)2 ,i + S 2 ,i,i) + C(3) (Si + 3 5 2 ) 



85 2 



+ 



24 So 



(n + 1) 3 (n + 1) 



+ 



16 5, 



16 S u 



6 5i5 3 2 5 2 



n + 1) 2 (n + 1) 



+ 



, | 45 4 | 165i, 3 



n + 1 n + 1 n + 1 n + 1 



+ ^-^ + ^-65i5 4 -5 2 5 3 -75 5 + 105i, 4 + 452,3 
n+1 n+1 n+1 

+ 5*1,1,3 — 4 5i,3,i + 3 53,1,1 — 5i,2,2 + 2 52,1,2 — 52,2,1 — 5i, 1,2,1 + 5i,2,l,l 
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3.4 H(l 



4; 2, 3) and H(l 



3; 2, 4) 



The H integral enjoys the same type of flip symmetry as the ladder graphs and E. In order to 
analyse H(l, 4; 2, 3) we may thus start from H(2, 3; 1, 4) instead. We send point 4 to infinity and 
identify x\ = 0, x 2 = pi, x 3 = p 2 ,x$ = ki, xq = k 2 , x-j = k 3 . In momentum space notation: 



If k\,k 2 ~ P2', k 3 ~ px or ki,k 2 ~ p\\ k 3 ~ p 2 we find no-scale integrals. The contributing 
regions are thus the top and bottom problems — both of topology FA — and the two regions 
R 2 : k\, k 3 ~ p 2 ; ki ~ p\ (and its mirror image with k\ k 2 ), and R3 : fci, k 3 ~ p\\ k 2 ~ p 2 
(and the same with k\ <B- k 2 ). 

The leading 1/e 3 pole must be universal to all regions for the recombination of the logarithms 
into powers of log(w) to happen. We can use the trivial cases R 2 , R 3 to try and understand the 
structure of a fit. Once again, the two types of numerator terms A = kf + k 2 and B = —2(ki.k 2 ) 
lead to different structures: The A terms give a series like Y n ~ x S(n)/n m , the B term causes 
Y n ~ x S \n) j \n + l) m . The Mincer system is able to deal with the top problem up to and including 
y 19 , thus we obtain 20 equations on the coefficients in an ansatz for each sector. It turns out 
that in the logarithm terms and the ((3) bit of the A part only Y n ~ x S^nj/n 1 occurs, see below. 
Imposing this for the rational terms at u° as well we can comfortably fix all coefficients. 

With some hindsight (Si t i 7 i,i, and high powers of Si do not occur) the ansatz for 

the B part can be limited to 28 coefficients, for which there are 20 equations, so what can be 
done? Exchanging x\ <H- x 2 is an active transformation mapping our case to H(l, 3; 2, 4). The 
usual transformation ([7j) on the cross ratios must be followed by dividing out (1 — Y) 2 due to the 
higher weight of the integral. This leads to different behaviour under the map: We find a single 
series Y n ~ l S(n)/(n + l) m , where for the first time the m = cases also occur. The resulting 
ansatz for the rational terms is fairly large because there are many independent harmonic sums 
(or products thereof) with positive indices adding up to total transcendentally weight six. 

Hence a priori we cannot expect any constraint from point exchange. As a matter of fact the 
simultaneous existence of the two expansions does impose a few conditions: We do not find any 
condition on the coefficients for the B fit for the log 3 (u) terms of H(l, 4; 2, 3) = H(2, 3; 1, 4), but 
one at log 2 (w), two at log(u) and three in the rational part. Hence we can restrict our ansatz for 
the rational terms of the B series to five unknown constants. Supplemented by the guess that 
all coefficients are integer multiples of eight it was not hard to play through the possibilities in a 
likely range; again we find one solution, which is presented below. This " diophantine" problem 
may seem a weak constraint, yet the difference in complexity is absolutely striking between 
the preferrred solution which we display and any other random try (where one puts integer or 
half-integer guesses for the independent parameters and inspects the values of the dependent 
parameters) . 



H(pi,p 2 ) = 



d D hd D k 2 d D k 3 (kj - 2(h.k 2 ) + fcf) 



(48) 



kfkikKh - h)\k 2 - h^ih-p^ih-ptfih-p^ih - P2 f 





l,4;2,3)-> 



(49) 



n=l 



\og{u) (8 S 2 2 + 16 S 1 S 3 + 4 S 4 - 8 5 1>3 + 8 5 1>lj2 - 8 5 2>ljl ) + C(3) (-16 S 2 ) 



+ 24 S 2 S 3 + 32 + 16 S 5 - 32 S 1A + 16 Si, li3 - 16 S 3 ,i,i + 16 S h2>2 - 16 5 2j2 ,i + 
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00 yn-1 

Y — 

^— ' 71+ 1 



n=l 



-g log 3 («) 



+ -io g 2 H 



- log(u) 



8£ X 



4£i 
n + 1 

4£ 2 



(n + 1) 2 n + 1 



-2S( + 4S 2 

8S 1 S 2 + 12 £ 3 



8 5? 



24 £ 2 8£ 3 16 51,2 



(n + 1) 2 (n + 1) 2 n + 1 n + 
- 8 5 lil)2 + 8 S 2)lil ) + C(3) (J^j + 16 Sf - 32 £ 



^ - 8 - 16 £i £ 3 + 20 £ 4 + 8 £i, 3 



64 gi 32 £ 2 64 £ 2 32 £i £ 2 16 £ 3 64 S* 1|2 16 £f 24 £ 4 

+ (n + l) 4 + (n + l) 3 ~ (n + l) 3 + (n + l) 2 ~ (n + l) 2 ~ (n + l) 2 ~ (n + 1) + n + 1 

_32£^ _ 165^ + ]^i_ 24S2S3 _ 32SiS4 + 24S5 + 32 5 i 65 ^3 + 32 53,1,1 
n+1 n+1 n+1 



- 16 5i, 2 , 2 - 16 5 2 ,i, 2 + 16 5 2 , 2 ,i 
4,if(l,3;2,4)-> 

oo 



(50) 



n=l 



~ log 3 («) 



4 5i 



+ 2 lo s 2 H 



- log(u) 



■Si 



n + 1) 

4S 2 | 4SiS 2 
(n + 1) 3 (n + 1) 2 n + 1 
4S 2 



2 + -^r + -^T - 2S 1 S 2 + 2S 3 + 2S h2 ) 
1 n + 1 n + 1 / 

— 4 Si S 3 + 12 S 4 + 2 5i,2,i — 2 £2,1,1 



12 5*2 8 5*i S 2 8 5*3 



8SV 



+ 



2S 2 2 14SiS 3 



(n + 1) 3 (n + 1) 3 (n + 1) 2 (n + 1) 2 (n + 1) 2 n + 1 n + 1 
12>Si,3 2 5^1,2 2 52,1,1 



+ 



+ 4 S 2 S 3 - 2 £1 S 4 + 42 S 5 - 14 5, 4 - 10 S 2 ■ 



n+1 n+1 n+1 
+ 4 £1,1,3 + 4 5*1,3,1 — 8 £3,1,1 — 2 £2,1,2 + 2 £2,2,1 — 2 £1,1,2,1 + 2 £1,2,1,1 
32 £1 8 S 2 8 £ 



+ C(3) 



[n + 1) 2 n+1 n + 1 



+ 16 £1 £ 2 - 16 £ 3 - 16 £1 2 



+ 



64 £1 16 £f 



16 £ 2 16 £1 £ 2 32 £ 3 32 £1,2 44 £1 £ 3 12 £ 4 



(n + 1) 5 (n + 1) 4 (n + 1) 4 (n + 1) 3 (n + 1) 3 (n + 1) 3 (n + 1) 2 (n + If 
56 £1,3 + 4 5i,i, 2 _ 4 5*2,1,1 + 16 5 2 5 3 + 44 5i 5 4 + 20 5 5 52 5 M 28 5 2 , 3 



+ 



(n + 1) 2 (n+1) 2 (n + 1) 2 n+1 n + 1 n+1 n+1 n+1 
8 £1,1,3 8 £ 3 ,i,! 4 £2,1,2 , 4 £2,2,1 , 2 o2o2 4 o2 , ooo , 20 



n+1 n + 1 n+1 



H t^t 1 + 77 £1 £2 — 77 £3 + 22 £2 £ 4 + — £1 £5 

n + 1 3 3 3 



^2 10 62 

+ 108 £ 6 - 48 £i, 5 - 56 £ 2 , 4 + — £i, 1>4 + — £1,4,1 - — £4,1,1 

26 4 26 38 28 14 

+ y ^1,2,3 — 7j ^1,3,2 + y ^2,1,3 + y £2,3,1 — ^ £3,1,2 + y £3,2,1 

8 2 14 2 

— 4 £1,1,1,3 — 4 £1,1,3,1 + 8 £1,3,1,1 — 7j £1,1,2,2 — 7j £1,2,1,2 — g- £1,2,2,1 — ^£2,1,1,2 

- y £2,1,2,1 + 7j £2,2,1,1 + 2 £1,1,1,2,1 - 2 £1,1,2,1,1 
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4 CPWA analysis 

At tree level, the correlation function is given by products of free scalar propagators: 



z' 7\r2 i\2 „,2 \ 2 / „,2 ,,2 \ 2 /„,2 „,2 \ 2 1 

0^(1,2,3,4) ' * 



4(4tt 



2U 



MM] + MM + (MM 

x 12 x 34/ \ x 13 x 24/ \ x 41 x 23 



(51) 



( /1'7t2^4 \ 2 ™ 2 ™ 2 ™ 2 /-v»2 , > . 2 , >-» 2 ™ 2 , >-» 2 ™ 2 , >-» 2 ™ 2 

^/l ; \ x 12 x 23 x 34 x 41 x 12 x 24 x 34 x 13 x 13 x 23 x 24 x 41 



The loop-corrections to G 4 take a factorised form [TJ [2] : 

(1,2,3,4) = 2( / 5~ 1) x #(1,2,3,4) x F«(x,) (for / > 1) (52) 



(4tt 2 ) 4 

Here i?(l,2,3,4) is a universal, /—independent rational function of the space-time, Xj, and har- 
monic, yj, coordinates at the four external points 1,2,3,4: 

Dfl O Q A^i _ VuVhyuy 14 ( 2 2 _ 2 2 _ „2 2 \ 
-TL^±, z,, o, — 2 2 2 2 V x 13 x 24 x 12 x 34 x 14 x 23/ 
x 12 x 23 x 34 x 14 

2/l22/l32/242/34 / 2 2 22 2 2 \ 

+ ™2 „2 ™2 „2 l X 14 x 23 x 12 x 34 X 13 X 24J 
X 12 X 13 X 24 X 34 

2/132/142/232/24 / 2 2 22 2 2 \ 

+ ™2 „2 ™2 „2 l X 12 x 34 x 14 x 23 X 13 X 24J 
X 13 X 14 X 23 X 24 
4 4 4 4 4 4 

1 2/122/34 1 2/132/24 1 2/142/23 /co\ 
~ l ~2 2~ l ~2 2~'~2 2> V 0,3 / 
x 12 x 34 x 13 x 24 x 14 x 23 

while F^'(xi) are functions of X{ only, which are explicitly stated up to three loops in terms of 
the box integrals and E and H in formulae @, fl3]) in the introduction. 

In the OPE limit x 2 — > X\, x± — > x 3 the weight 2 operators at points 1,2 and 3,4 respectively 
fuse into an expansion in terms of operators "exchanged" between the two halves of the four- 
point function. The exchanged operators carry twist (dilatation weight - spin), spin and 577(4) 
quantum numbers. Since we are fusing two operators in the 20' representation, the exchanged 
operators must carry one of the representations in the product 20' x 20' = 1 + 15 + 20' + 84 + 
105 + 175. The correlator as written in the last two formulae has six "channels" distinguished 
by the y variables pertaining to the internal symmetry group. It has been worked out in [16] 
which linear combination of these channels correspond to the exchange of operators in a given 
representation. In particular, if we label the channels according to 

G 4 = 2/122/34 ^1+2/132/24 A 2 +2/142/23 A z +yl 2 y34yi3y2A A ^+yi2y34yuy23 A 5+yi3yuyuyL A 6 (54) 

then the 20' exchange corresponds to A 2 + A 3 + 1 At + 1 A 5 + 1 A 6 . Both at tree- and at loop-level 
we find that the leading power singularity is l/(x 2 2 x| 4 ) coming from A4,A 5 : 

^Gf = ot J^f~ 2 \ 2 ^ 2 * (55) 



zi2,z34->o 3 (An 2 ) A x\ 2 x 3A x\ 3 x 



2 
21 



Y 



+ 1 ) +2 (r 1 ^- 1 - Y )f^ lim n «'^W|+o(«) 



1=1 
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where a is the effective coupling. Powers of u can be discarded because they correspond to 
higher twist. On the other hand, the expansion in the Y variable is associated to the spin of the 
exchanged twist two operators. 

In terms of the elementary field the operators we discuss are schematically realised as = 
tiiWD^ 1 . . .D^W). The positioning of the trace-free symmetrised Yang-Mills covariantised 
derivatives = <9 M + % g [A 11 , •] on the two W fields decides about the anomalous dimension 
of the operators. All these composites have twist two because there are two scalar elementary 
fields of dimension one carrying the derivatives. The spin operator in this class is the protected 
primary O. At spin one we only find d^O, a "conformal descendent" of O. By tree-level 
orthogonalisation one sees that there is one new primary operator at every even spin, all other 
combinations of derivatives give descendents of operators with lower spin. The descendents do 
show in the OPE, but their occurrence is statically linked to that of the primary field. One 
can resum the contribution of the descendents into a "conformal partial wave" labelled by the 
primary operator, see [SI EH EZ] and references therein: 

cpwa(s) = N(s) v? /2 Y s 2 F 1 (s + 1 + 7/2, s + 1 + 7/2, 2 + 2 s + 7; r) (56) 

When expanding in the effective coupling a the cpwa furnish logarithms to be matched on those 
of the in fl56|) . Our task is thus to solve 

1- 2 ( — 1 - y] V lim a 1 F m (xi) | = V cpwa(s) (57) 

v 7 1=1 / s=0 



Y 



for 7(5), N(s). As it should, the normalisation of the cpwa turns out to be zero if the spin is odd. 
For even spin we find 




r (2 s + 1 + 7) 

* + a ci, 2 + a 2 (C(3) C2,i + c 2 , 4 ) + a 3 (((5) c 3 ,i + C(3) c 3 , 3 + c 3 , 6 ) + 



where 



g 2 N 



Air 



2 ■ 



7 2 + (S 1 (2 S )-5 1 ( S )) 7 3 + ... , 6 3 = 7 3 + -- - • (59) 



The 6j cancel the explicit dependence of the T ratio on ( values; the first factor in (1581) is meant 
to be purely rational. The anomalous dimension depends on the spin and has an expansion 
7 = a 71 + a 2 72 + a 3 7 3 + . . . 

71 = 251, (60) 

72 = -5 , _ 3 -2 5 , _ 2 5 , 1 -25 , i5 , 2 -5 , 3 + 2 5 , _ 2jl , (61) 

73 = 3 5_ 5 + 8 S_ 4 ^ + S? 2 S 1 + 6S- 3 Sl + S^S 2 + AS^ 2 S 1 S 2 + 2S 1 S 2 (62) 

+ 2 5_ 2 S 3 + 2 5? S 3 + 5 2 5 3 + 3 Si S 4 + S 5 - 6 5-4,1 - 12 Si S-3,1 - 6 S_ 3>2 
- 4 5 X 2 5_ 2 ,i - 2 5 2 S_ 2 ,i - 10 Si S_ 2 ,2 - 6 5-2,3 + 12 5-3,i,i + 16 Si 5_ 2 ,i,i 
+ 12 5_2,i,2 + 12 5-2,2,1 — 24 5_2,i,i,i • 
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In these formulae all harmonic sums depend on the argument s. Our result is in complete 
agreement with the literature [10] . 

The tree-level normalisation is easily recognised to be the ratio 2 (s!) 2 /(2s)! which we take 
out of the entire normalisation factor. A fit of the one-loop normalisation on harmonic sums is 
then in fact possible, but only if in addition S(2s) are taken into account. The result for the 
coefficients is compatible with promoting the the factorials at tree to the T functions in ( 158]) ; 
this mimics the Y dependent part of the cpwa. At two and three loops the T functions correctly 
incorporate all S(2s) terms. Below we state our results for the coefficients in the second line of 
( I58p . In these formulae all harmonic sums have argument s once again. 

ci, 2 = -5*2, (63) 

C2,i = 3S 1; (64) 

5 5 
c 2,4 = 2 5-4 + S_ 2 + 2 5_3 Si + >S_2 5*2 + S 2 + 2 Si S3 + — S4 — 2 S_3i — S_ 2 ,2 ~~ 2 S13 , (65) 

C3,i = _ y5i, (66) 

c 3 , 3 = -3S_ 3 -10S^Si + ^S!-QSiS 2 -^S 3 + QS^ u (67) 

5 41 222 

c^fi = — 11 S_6 + — S_ 3 — 5 »S_4 S^2 — 2~ 5-5 Si — S-3 S-2 Si — 5 S_4 Si — 2 S_ 2 Si (68) 

1 c c3 c c o2c inc cc oc c2 c3 c c 

3 — 2~ 2 _ _ _ — 3~ 

15 13 

— 8 S_2 Si S3 — 6 Si 52 53 — 3 S3 — 3 S*_2 54 + 9 Si S4 — 4 S 2 54 H — — 5i S5 — — 5g 

+ 14 S_ 5 ,i + 11 5i S_ 4 ,i + 9 5_ 4 ,2 - 12 5i S_ 3 ,- 2 + 10 5_ 2 5-3,1 - 4 Si S_ 3 ,i 

+8 S 2 5_ 3 ,i + 4 5 X S_ 3 , 2 + 9 5_ 3 ,3 - 10 5_ 3 5_ 2>1 + 14 S_ 2 5i S_ 2 ,i - ? Sj 5 S_ 2 ,i 
20 

+ ASiS 2 S-2,1 + —S 3 5_ 2 ,i + 10 S\i + 10 S_ 2 5_ 2 ,2 - 6 5 X 2 S_ 2 , 2 + 6 5 2 S_ 2 , 2 
o 

+ 6 5i S_ 2 , 3 + 11 5_ 2 ,4 - 6 5 2 5i, 3 -A Si S M - 4 S 1>5 + 4S a 5 2 , 3 + 4 S 2 , 4 - 12 S_ 4 ,i,i 
+ 8 5*i 5-3^,1 — 2 S_3,i )2 — 2 5-3,2,1 — 24 Si S_ 2 ,_ 2 ,i — 20 S_2 5-2, 1,1 + 16 S 2 5_2,i,i 

— 8 S 2 S— 2,1,1 + 16 Si S_ 2 ,l,2 — 6 S_2,l,3 + 16 Si S_ 2 , 2 ,l + 4 S_ 2 ,2,2 — 6 S_2,3,l — 4 Si Sl,l,3 

— 8 Si,i,4 + 8 Si^^ — 8 S-3,1,1,1 — 48 Si S_ 2 , 1,1,1 — 20 S_2, 1,1,2 — 20 S_2, 1,2,1 — 20 S_ 2 , 2) i,i 
+16Si i i i i i3 + 64S_2 ,1,1,1,1 



5 Conclusions 

In a double concidence limit x 2 — > x±, £4 — > X3, the four-point function of stress tensor multiplets 
T is reduced to an OPE (operator product expansion) T{x\)T{x2) = ^2 s c(s,x) C (^i)) an d 
similarly at the other end, so that one obtains (the sum over) the square of the structure con- 
stants c with the two-point function of the exchanged operator 0^ s \ This is a scheme invariant 
combination from which one can read off the structure constants if the two-point function in the 
middle is assumed to be normalised to one. 
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Not only the primary operators but also their conformal descendents (x-derivatives) are ex- 
changed. The descendents are usually put together with the primary fields in conformal blocks 
called "conformal partial waves" (cpwa). We have derived explicit results for the twist two 
operators in the 20' representation of S77(4): Their anomalous dimensions come out as linear 
combinations of harmonic sums in full agreement with the literature [9j[T0]. What is more, also 
the constants multiplying the twist two cpwa are elaborated in terms of harmonic sums. 

We have stopped short of predicting the structure constants from these results because the 
absolute normalisation of the cpwa is not knowrJl. Our results for the constant terms naturally 
factors into two pieces. It is tempting to associate the first of these factors with the normalisation 
of the cpwa and the second with the structure constants. 

The fact that the entire result is expressed in terms of harmonic sums is a clear hint at an 
integrable systems explanation, c.f. [TH [12]. This issue will be addressed in future work; we are 
confident that the normalisation question will be understood if such an interpretation is found. 

The asymptotic expansion of the individual conformal integrals is given in terms of harmonic 
sums with positive indices only, and products thereof with negative powers of their argument. 
These results should help with the construction of an explicit expression for the unknown integrals 
in terms of special functions of the polylogarithm type. Interestingly, the cpwa decomposition 
leads to formulae in terms of harmonic sums only, but here the sums can have negative indices. 
S-x does not occur, and in most of the higher sums only the outermost index can be negative 
(two exceptions). 

In deriving the asymptotic series by expansion by regions we have met a number of structural 
properties, i.e. that given numerator terms lead to expansions that can be matched on the 
distinct structures Y n ~ l S(n) / n m , y n_1 S(n) /(n + l) m with or without the m = case. We have 
used conformal symmetry to make four-point integrals into three-point integrals. The latter are 
generic by inspection; one may wonder whether any three-point integral can be written as a spin 
expansion in terms of harmonic sums, or whether the examples here are somehow specific to the 
J\f = 4 SYM theory. 

Obviously, our work can be extended to the twist three, four, ... trajectories corresponding to 
powers of the second cross ratio u. We expect that the coefficients will pick up a second parameter; 
it remains open for the moment whether the rational number in front of each harmonic sum will 
simply start to depend on the twist, whether each trajectory is completely different, or indeed if 
the Euler-Zagier sums are not sufficient to express the complete expansion. 

Last, on very many occasions — but not always — the coefficients of a set of harmonic sums 
related by index permutations add up to zero in our formulae. This hints at the existence of a 
special basis w.r.t. which the results would take a much simpler form. 
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6 Appendix: The 20 channel in terms of harmonic sums 



In the four-point correlation function at two and three loops the integrals come with very specific 
rational factors depending on Y . It is an interesting question how the expansions in terms 
of harmonic sums over powers of their arguments change by including the Y expansion of the 
rational factors, and which type of expansion will lead to a similar result if multiplied by such 
factors. In the 20' channel we obtain a fit of the same type as for the asymptotic expansions of 
the individual integrals even if the factor \/{\ — Y) — \— Y from the R polynomial projected onto 
the 20' representation is included. We display the result as an illustration of the universality of 
the basis of harmonic sums we were using, but also in the hope that the formulae may be useful 
in an attempt on deriving an explicit result for the correlation function in terms of harmonic 
polylogarithms and related functions. The tree-level 1/(1 — Y) + 1 + 0(u) is thus followed by 
the loop correction 2 Y^iti a 1 + 0(u) with the twist 2 contributions 



where 



71=2 



log(u) (-— + -S 1 )+--S: 



2n 



/( 2) = Y jY n 



n=2 



\0g\u) C 2 , 2 + log(u) C 2 , 3 + C 2 ,4 + C(3) #2,1 



(69) 
(70) 



and 

/ (3) = 

with 



a 



2,2 



a 



2,3 



G 



2,1 



- + 351, 
n 



7> — 9 + 7 ^1 + 7 ^2 
2n 2 4 4 



n° 



S\ 5*2 3 3 
2^ + 2^ 2 ^ 2 2 3 ' 



6 Si 

7 + 4 + 



n 



rr 



2n 2 
a 



+ 



S\ Si 



2SV 



n 2n n 

I 13 5 

S\ Sip + - S 2 + 3 Si S3 + — S*4 — — Si^ + 3 Si, 1,2 



(71) 
(72) 
(73) 
(74) 



n=2 



log 3 («) C 3 , 3 + log 2 ( M ) C 3;4 + log(u) {C 3 , 5 + C(3) L> 3 ; 2 ) + C 3;6 + C(3) D 3 . 3 + C(5) D 3; 

(75) 



25 



3:1 



- 25^, 



^3;2 = 



— + 3 Sf + 3 S 2 , 



14 10 5i 2 Sj 2S 2 2 3 14 

-^3;3 — — ? 5 1 1" >->l — 8 S 2 Si — D 3 + 8 612 , 

n 6 n z n n 3 3 



(76) 
(77) 
(78) 
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3:4 



7>' J 



C C2 
~~ 3 



n 



An 2 An 2 



S2 Si 5*2 
2n 



^3 
2n 



2n 



o2 o C C C C 

D 1 D 2 — Oi D 3 — - D 2 — - 04 — - Di i3 , 



30) 



a 



3;5 



30 6 Si Sf 



72 1 



2n 3 2n 
3 Sj^ 5*2 Si 53 3 S2 



5 S2 , 5 S\ S2 5 S3 



2n 2 



2n 2 



2n 
9 



An 
13 

~4~~ l ~ u ■ 4 

+ 5*1,2,2 + 6 Si 1 1,2 



H — 7- S 2 S3 + — Si S 2 



An 

> 29 



S4 2 5*1,3 1 C 2 C 

0^1 ^l,: 

n 71 2 



2 5*1,2 3 5i 5i,2 
n 2 2n 

Si 5i,3 — 2 5*2 5*1,2 



17 

— 5i 54 + — 52 53 



6 5*5 — 5 5*i 4 + S2 ; 



Si, 



1,3 



51) 



3:6 



60 12 Si 6 5 : 



72 " 



3 5*. 



n° 

2 Sf S2 



2n s 
5i 2 5 3 



+ 



2n 2 



1 sf 


3 5*i 5*2 


65*3 


n 4 


n 3 


n 3 


6 5i5 3 


8 5*4 


10S 1)3 


n 2 




n 2 



+ 



4 Si 2 Si 5 



1.2 



6 S2 5i,2 2 Si Si, 3 



5 Si S^ 



13 Si S4 3 S5 16 Si, 4 2 S2,3 



72 



2n 



72 



72. 



72 



72 

3 Si, 2, 2 

n 



72 

125 
IT 



2n 



13 



Sq + TTT So 



12 



21 1 2 

+ 3 Si S2 Si,2 + — S2 Si,3 + 2 S3 Si, 2 

1 



11 



\2) 



2 5*i Si,3 + 3 Si Si,4 — Si S2,3 — — S2 Si, 1,2 



41 



~~ 2 ^i 5i,i,2 — Si Si, 1,3 — 3 Si Si, 2,2 + 4 Si Si, 1,1,2 — - S x S 2 — 9 Si S2 S3 — — S2 S4 
7 17 

— — S 3 — — S 1 S4 — 19 Si S5 + S 1 2 + 24 Si, 5 + 6 S2,4 + 2 Si, 1,4 — Si, 2, 3 



6 Si,3,2 — 3 Si, 1,1,3 + Si, 1,2,2 — 10 Si, 1,1, 1,2 



References 

[1] B. Eden, A. C. Petkou, C. Schubert and E. Sokatchev, "Partial nonrenormalization of the 
stress tensor four point function in M = A SYM and AdS / CFT," Nucl. Phys. B 607 
(2001) 191, hep-th/0009106[ 

[2] B. Eden, P. Heslop, G. P. Korchemsky and E. Sokatchev, "Hidden symmetry of four- 
point correlation functions and amplitudes in N=4 SYM," Nucl. Phys. B 862 (2012) 193, 
I arXiv: 1108.35571 [hep-th] . 

[3] B. Eden, P. S. Howe, C. Schubert, E. Sokatchev and P. West, "Four point functions in 
Af = A super symmetric Yang-Mills theory at two loops, " Nucl. Phys. B557 (1999) 355-379, 
hep-th/98 11172] "Simplifications of four point functions in N ' — A supersymmetric Yang- 
Mills theory at two loops," Phys. Lett. B466 (1999) 20-26, ^^1:11/9906051) F. Gonzalez- 
Rey, I. Y. Park and K. Schalm, "A Note on four point functions of conformal operators in 
Af = A superYang- Mills, " Phys. Lett. B448 (1999) 37-40, |hep-th/98lTT55| 



21 



[4] B. Eden, C. Schubert and E. Sokatchev, "Three loop four point correlator in M = 4 SYM, " 
Phys. Lett. B482 (2000) 309-314, |hep-th/0003096[ 

[5] M. Bianchi, S. Kovacs, G. Rossi and Y. Stanev, "Anomalous dimensions in M = 4 SYM 
theory at order g**4," Nucl. Phys. B584 (2000) 216-232, [hep-th/0003203| 

[6] L. F. Alday, B. Eden, G. P. Korchemsky, J. Maldacena and E. Sokatchev, "From corre- 
lation functions to Wilson loops," JHEP 1109 (2011) 123. larXiv: 1007.3 243 [hep-th] ; 
B. Eden, G. P. Korchemsky and E. Sokatchev, "From correlation functions to scattering 
amplitudes", JHEP 1112 (2011) 002. larXiv: 1007.32461 [hep-th]; "More on the duality 
correlators /amplitudes, Phys. Lett. B 709 (2012) 247, larXiv: 1009T2 488 [hep-th] . 

[7] B. Eden, P. Heslop, G. P. Korchemsky and E. Sokatchev, "Constructing the correlation 
function of four stress-tensor multiplets and the four-particle amplitude in Af = 4 SYM," 
Nucl. Phys. B 862 (2012) 450, arXiv: 1201 .~5529~ [hep-th] . 

[8] I. T. Todorov, M. C. Mintchev and V. B. Petkova, "Conformal Invariance in Quantum 
Field Theory," Pisa, Italy: Sc. Norm. Sup. (1978) 273p. 

[9] F. A. Dolan and H. Osborn, "Superconformal symmetry, correlation functions and the 
operator product expansion, " Nucl. Phys. B 629 (2002) 3, |hep-th/0112251[ 

[10] A. V. Kotikov, L. N. Lipatov and V. N. Velizhanin, "Anomalous dimensions of Wilson oper- 
ators in N=4 SYM theory, " Phys. Lett. B 557 (2003) 114, |hep-ph/0301021 ; A. V. Kotikov, 
L. N. Lipatov, A. I. Onishchenko, V. N. Velizhanin, "Three loop universal anomalous di- 
mension of the Wilson operators in N=4 SUSY Yang-Mills model," Phys. Lett. B595 
(2004) 521-529, hep-th/0404092| 

[11] M. Staudacher, "The Factorized S-matrix of CFT/AdS," JHEP 0505 (2005) 054, 
hep-th/041 2188[ N. Beisert and M. Staudacher, "Long-range psu(2,2 — 4) Bethe Ansatze 
for gauge theory and strings, " Nucl. Phys. B 727 (2005) , |hep-th/0504190[ B. I. Zwiebel, 
"M = 4 SYM to two loops: Compact expressions for the non-compact symmetry algebra of 
the sm(1,1|2) sector," JHEP 0602 (2006) 055, lhep-th/0511109| 

[12] J. Escobedo, N. Gromov, A. Sever and P. Vieira, "Tailoring Three-Point Functions and 
Integrability," JHEP 1109 (2011) 028, larXiv: 1012 . 24751 [hep-th]; N. Gromov and 
P. Vieira, "Quantum Integrability for Three-Point Functions. " larXiv: 1202.4 103 [hep-th]. 

[13] V. A. Smirnov, "Applied asymptotic expansions in momenta and masses, " Springer Tracts 
Mod. Phys. 177 (2002) 1. 

[14] K. G. Chetyrkin and F. V. Tkachov, "Integration by parts: The algorithm to calculate beta 
functions in 4 loops," Nucl. Phys. B 192 (1981) 159; S. G. Gorishny, S.A. Larin, L.R. Sur- 
guladze and F.V. Tkachov, "Mincer: Program for multiloop calculations in quantum field 
theory for the schoonschip system," Comput. Phys. Commun. 55 (1989) 381; S. A. Larin, 
F. V. Tkachov and J. A. M. Vermaseren, "The FORM version of MINCER, " Preprint 
NIKHEF-H/91-18 (Amsterdam 1991). 



22 



[15] N. I. Usyukina and A. I. Davydychev, "An Approach to the evaluation of three and four 
point ladder diagrams," Phys. Lett. B 298 (1993) 363; "Exact results for three and four 
point ladder diagrams with an arbitrary number of rungs," Phys. Lett. B 305 (1993) 136. 

[16] G. Arutyunov, S. Frolov and A. C. Petkou, "Operator product expansion of the lowest weight 
CPOs in N=4 SYM(4) at strong coupling," Nucl. Phys. B 586 (2000) 547 [Erratum-ibid. 
B 609 (2001) 539], hep-th/0005182| 

[17] F. A. Dolan and H. Osborn, "Conformal four point functions and the operator product 
expansion," Nucl. Phys. B 599 (2001) 459, |hep-th/0011040| 



23 



